Quantum disorder and Griffiths singularities in bond-diluted two-dimensional 

Heisenberg antiferromagnets 
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We investigate quantum phase transitions in the spin- 1/2 Heisenberg antiferromagnet on square 
lattices with inhomogeneous bond dilution. It is shown that quantum fluctuations can be contin- 
uously tuned by inhomogeneous bond dilution, eventually leading to the destruction of long-range 
magnetic order on the percolating cluster. Two multicritical points are identified at which the 
magnetic transition separates from the percolation transition, introducing a novel quantum phase 
transition. Beyond these multicritical points a quantum-disordered phase appears, characterized 
by an infinite percolating cluster with short ranged antiferromagnetic order. In this phase, the 
low-temperature uniform susceptibility diverges algebraically with non-universal exponents. This is 
a signature that the novel quantum-disordered phase is a quantum Griffiths phase, as also directly 
confirmed by the statistical distribution of local gaps. This study thus presents evidence of a genuine 
quantum Griffiths phenomenon in a two-dimensional Heisenberg antiferromagnet. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Mg, 64.60.Ak 



I. INTRODUCTION 

Quantum phase transitions in low-dimensional quan- 
tum Heisenberg antiferromagnets (QHAFs) have been 
the subject of extensive investigations during the last 
two decadesii. The transition from a renormalized clas- 
sical to a quantum-disordered state can be triggered by 
various parameters, including lattice dimerization, frus- 
tration and appUed field. More recently, special inter- 
est has focused on phase transitions driven by geometric 
randomness of the lattice^, including site and bond dis- 
order. Strong geometric disorder not only breaks trans- 
lational invariance and perturbs the ground state of the 
pure system, but it can also destabilize renormalized clas- 
sical phases with long-range order (LRO) and drive the 
system to novel disordered phases. 

Various one-dimensional QHAFs with bond disorder 
have been found to display unconventional quantum 
phases^. For example, the undimerized QHAF chain 
is driven into a random- singlet phased with algebraically 
decaying spin-spin correlations. The low temperature 
susceptibility diverges as l/[Tlog^T] in this phase, inde- 
pendent of the details of the bond disorder. In contrast, 
a dimerized chain shows a quantum Griffiths phased be- 
yond a critical disorder strength, with exponentially de- 
caying spin-spin correlations and a non-universal power- 
law divergent susceptibility. In two dimensions, the clean 
QHAF develops antiferromagnetic (Neel) LRO at zero 
temperature^. This introduces the intriguing possibility 
of a genuine order-disorder transition, driven by lattice 
randomness, i.e. from the Neel phase into one of the 
above unconventional disordered phases^. 

Site and bond dilution of the square lattice QHAF 
have been the focus of several recent studies^i^iiSii^ii^, 
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motivated by experiments on antiferromagnetic cuprates 
doped with nonmagnetic impurities^iiSiii. From a geo- 
metric point of view, bond and site dilution reduce the 
connectivity of the lattice, ultimately leading to a per- 
colative phase transitioniS beyond which the system is 
broken up into finite clusters. In a classical spin sys- 
tem, this percolation transition is coupled to a magnetic 
transition with the same critical exponents since sponta- 
neous magnetic order cannot survive beyond the perco- 
lation threshold. In a quantum spin system, on the other 
hand, a progressive reduction of the lattice connectivity 
enhances quantum fluctuations in a continuous fashion, 
raising the possibility of quantum destruction of mag- 
netic order before the percolation threshold is reached. 
However, recent studies of homogeneously site- and bond- 
diluted QHAFs on the square lattico^i^d^i^° found that 
the magnetic transition takes place exactly at the perco- 
lation threshold, due to the fact that the percolating clus- 
ter at threshold shows LRO''. The critical exponents of 
the correlation length and of the order parameter, v and 
/3, are found to take their classical percolation value3^*2i. 
Therefore, in this case the magnetic transition is com- 
pletely dominated by classical percolation^'^. 

Alternative to the above classical percolation picture 
is the quantum percolation mechanism, recently demon- 
strated in a number of model system322i2iiSS42&. This sce- 
nario is based on the fact that spins involved in locally 
strongly fluctuating quantum states, such as dimer sin- 
glets and resonating valence bonds (RVBs), are weakly 
correlated with the remainder of the system. In a random 
network of spins, the local strongly fluctuating states cre- 
ate weak links with small spin-spin correlations. If these 
weak links are part of the backbone of the percolating 
cluster, they can prevent the percolating cluster from de- 
veloping long-range order. Therefore, if lattice dilution 
favors the local formation of such states, it is possible to 
drive the system towards a quantum disordered state be- 
fore the percolation threshold is reached, thus decoupling 
percolation from magnetic ordering. 
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In this work we numerically scrutinize this quantum 
percolation scenario in the S — 1/2 QHAF on a square 
lattice with inhomogeneous bond dilution^i. The inho- 
mogeneous character of lattice randomness is the key in- 
gredient to enhance quantum fluctuations, favoring the 
local formation of RVB/dimer-singlet states. Based on 
large-scale Quantum Monte Carlo (QMC) simulations, 
we observe that the classical percolation scenario is pre- 
served for moderate inhomogeneity, up to a multicrit- 
ical point beyond which the magnetic transition devi- 
ates from the percolation threshold. At this point, the 
quantum percolation scenario sets in, leading to a novel 
quantum-disordered phase. Such a phase has very un- 
conventional features, i.e. a finite correlation length but 
a gapless excitation spectrum, and a diverging uniform 
susceptibility at zero temperature with a non-universal 
exponent. These signatures allow us to identify this 
phase with a genuine two-dimensional quantum Griffiths 
phaseW^. 

This paper is organized as follows. In Sec. ^Ithe model 
is described, and in Sec. IIIII some technical aspects of 
the simulations are reviewed. In Sec. IIVI we discuss in 
detail the fate of antiferromagnetic order on the perco- 
lating cluster upon tuning the inhomogeneity of bond 
disorder. The complete phase diagram of the model is 
discussed in Sec. |3 Sec. IVII deals with the emergence of 
the quantum-disordered phase in the limit of randomly 
coupled ladders, whereas Sec. IVIII is dedicated to the 
novel ground state properties of the quantum-disordered 
regime. Conclusions are drawn in Sec. IVlllI 



II. INHOMOGENEOUS BOND DILUTION OF 
THE QHAF ON THE SQUARE LATTICE 

We investigate the ground state and thermodynamics 
of the 5 = 1/2 QHAF on a two-dimensional square lattice 
with inhomogeneous bond dilution. The Hamiltonian of 
this system is given by 

^ = E 4''^ s,-s, + j J2 4'"^ ■ (1) 

{ij)eD (im)G-L 




FIG. 1: (color online) (a) Decomposition of the square lattice 
into dimer {D) and ladder (L) bonds; (b) randomly coupled 
dimer limit; (c) randomly coupled ladder limit. 




FIG. 2: (color online) Inhomogeneous bond dilution favors 
the formation of local strongly fluctuating quantum states. 
The green ellipses indicate such local dimer-singlet, plaquette- 
singlet and ladder-like RVB states: (a) Pl > Pd', (b) Pd > 
Pl. 



The sums run over the dimer (D) and ladder (inter- 
dimer) (L) bonds, as indicated in Fig. QJa). The ex- 
change couplings are taken to be equal and antiferro- 
magnetic (J > 0). ea{a = D,L) is a random variable 
drawn from a bimodal distribution taking values 1 ('on') 
or ('off'). The inhomogeneity of this model is of purely 
statistical nature, stemming from the different probabil- 
ities of assigning the 'on'-state {ci — 1) to the dimer 
bonds versus the ladder bonds. For D bonds this prob- 
ability is denoted by p{eD = 1) = Pd, whereas for L 
bonds it is denoted by p{eL ~ 1) = Pl- The special 
case Pd — Pl represents the previously studied homoge- 
neous bond-dilution problem^'l® . In the general inhomo- 
geneous case Pd ^ Pl, two limits are noteworthy (see 
Fig. nb)-(c)): the case Pd = 1, Pl G [0, 1] corresponds 



to a system of randomly coupled dimers, and the oppo- 
site case Pl = I, Pd G [0, 1] corresponds to randomly 
coupled ladders. The degree of inhomogeneity can be 
parametrized by introducing a variable 

e = arctan (j—^^ (2) 

such that the limiting cases of randomly coupled lad- 
ders, homogeneous bond dilution and randomly coupled 
dimers correspond to 6' = 0, 7r/4, and tt/2 respectively. 

It is evident that the inhomogeneous nature of bond 
dilution in the lattice enhances quantum fluctuations in 
the magnetic Hamiltonian. In the strongly inhomoge- 
neous limits 6^0, tt/2 the structure of the perco- 
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lating cluster is geometrically built from weakly coor- 
dinated segments of ladders (0 — > 0), and weakly coordi- 
nated dimers {6 — > 7r/2), e.g., finite necklace-like struc- 
tures and finite randomly decorated chains (see Fig. |2J). 
Both, antiferromagnetic Heisenberg ladders'^*' and dec- 
orated chains^ii^Si^, have quantum- disordered ground 
states with a finite correlation length ^o- In antiferromag- 
netic ladders the ground state has a RVB nature'^*, and 
for the necklace it has a dimer-singlet nature^i*^. When 
segments of such structures become part of a larger clus- 
ter, as it is the case in the inhomogeneous percolation 
model, they locally retain ground state properties simi- 
lar to the thermodynamic limit if their length I is large 
compared to the correlation length ^o- This represents 
the core mechanism of non-trivial enhancement of quan- 
tum fluctuations through bond dilution, leading to novel 
quantum disordered phases. 



III. NUMERICAL METHODS 

The purely geometric problem of inhomogeneous per- 
colation is studied using a generalized version of a highly 
efhcient classical Monte Carlo algorithm.^^. In turn, the 
quantum magnetic Hamiltonian Eq. (|5J is investigated 
by Stochastic Series Expansion (SSE) QMC simulations 
based on the directed- loop algorithm'^'^. The ground state 
properties are systematically probed by efficiently cool- 
ing the system down to its physical T = behavior via a 
successive doubling of the inverse temperature /3 (Ref. 0). 
This approach is necessary since in two-dimensional (2D) 
diluted systems 2D correlations are often mediated by 
narrow quasi-one-dimensional links, such that the tem- 
perature scale at which these correlations set in is much 
lower than in the clean 2D case. 

The scaling analysis of the quantum simulation results 
is carried out in two complementary ways^. In the first 
approach, the lattice size is fixed at L x L, and the en- 
tire system is simulated, including percolating and non- 
percolating finite clusters. In the second approach, a fun- 
damental theoretical tool is used to ascertain the presence 
or absence of order in the system: we focus exclusively 
on the percolating cluster, i.e. discarding the finite clus- 
ters that cannot carry long-range order. On a finite Lx L 
lattice with periodic boundary conditions the percolating 
cluster is identified as the largest cluster, winding around 
at least one of the two spatial dimensions. The advantage 
of this approach is that a single run of QMC simulation 
evaluates the observables on both the full lattice and the 
largest (percolating) cluster. At the same time the largest 
cluster has non- negligible size fluctuations (of order L), 
leading to a signiflcant error in the QMC data. An alter- 
native approach, eliminating these size fluctuations, is to 
grow the percolating cluster freely from an initial seed up 
to a fixed size Nc but without any restriction of lattice 
boundaries. Starting from a single occupied bond, one 
visits the six neighboring bonds, activating them with 
probability Pjj (Pl) if they are D (L) bonds, etc. The 



procedure is repeated until no new bonds are generated. 
The cluster geometry is accepted only if its final number 
of sites matches exactly the required Nc. A drawback of 
this procedure is that the cluster build-up becomes time 
consuming for large Nc and away from the percolation 
threshold because of the large rejection rate. 

In the full-lattice studies, L values up to 64 were used. 
Within the /3-doubling scheme, inverse temperatures as 
high as f3 — 256L have proven to be necessary to observe 
the physical T = behavior. For the studies on fixed-iVc 
clusters we considered Nc values up to 2048 and inverse 
temperatures as high as /3 = 8Nc. For each lattice size 
and for each point P={Pd,Pl), 10^-10'^ realizations of 
the diluted lattice/percolating cluster were generated in- 
dependently to obtain disorder- averaged observables. 



IV. MAGNETISM ON THE PERCOLATING 
CLUSTER 

In this section, QMC results are discussed which ad- 
dress the evolution of antiferromagnetic order on perco- 
lating clusters upon tuning the nature of bond dilution 
from homogeneous {0 = 7r/4) to inhomogeneous {9 — > 0, 
7r/2). A fundamental observation^ is that the presence 
or absence of magnetic order in a diluted system and 
the nature of the transition from order to disorder are 
determined by the magnetic behavior of the percolating 
cluster. In two dimensions, only one percolating clus- 
ter can exist in the system. We denote with (Mc) the 
disorder- averaged ((...)) staggered magnetization per site 
of a percolating cluster, which is estimated by 




Here the summation ^ is restricted to sites contained 
in the percolating cluster, and Nc is the total number 
of these sites. This is to be contrasted with the overall 
magnetization of the diluted lattice with N — LP' sites. 




Given that only the percolating cluster contributes to 
long-range order in the thermodynamic limit, one finds 
that for N ^ oo 




Moreover, exploiting the self- averaging property of the 
size distribution of the percolating cluster—, close to the 
percolation threshold one can use 

(Af(P)) « A\P - MciP), (6) 
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FIG. 3: (color online) Upper panel: classical phase dia- 
gram for the inhomogeneous bond-percolation problem on the 
square lattice. The angle 9 parametrizing the critical curve is 
indicated. Lower panel: scaling plots for the network strength 
Nc/N in the two strongly inhomogeneous cases Pd =0.9 and 
Pl = 0.9. An excellent data collapse is realized, using the 
2D percolation exponents /3 = 5/36 and v — 4/3, and with 



0.4007 (Pd = 0.9) and P^'^'j = 0.0107 (Pl = 0.9) 



where P = {Pd,Pl) is the control parameter of the 
inhomogeneous percolation problem, and P — Pc"^^^ 
defines the classical critical line of percolation thresh- 
olds in the {Pd,Pl) space. In Eq. JB)) we have used 
the critical behavior of the so-called network strength 
{Nc/N) A\P - Pi'^'V, where A is a constant ampli- 
tude, and P = 5/36 is the critical exponent of classical 
percolationi^. 

From Eq. © it is evident that, if MdPc"'^) ^ 0, the 
dilution-driven transition from magnetic order to disor- 
der coincides with the classical percolation transition, 
both in terms of its location and of its critical exponents. 
Therefore only the absence of antiferromagnetic order on 
the percolating cluster can lead to a decoupling between 
the magnetic transition and the geometric percolation 
threshold. 

The critical curve P = Pc"^^^ , determined from classical 
Monte Carlo simulations, is shown in Fig. |21 Despite 
the inhomogeneous nature of percolation for P^ ^ -Pd, 
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FIG. 4: (color online) Phase diagram of the spin-1/2 QHAF 
on the inhomogeneously bond-diluted square lattice. 



excellent scaling of the classical MC results is observed 
using the 2D percolation exponents^^. This demonstrates 
that the universality class remains unchanged along the 
critical curve. 

Starting from the classical Monte Carlo result, the evo- 
lution of the order parameter can be tracked on the 
percolating cluster {Mc{Pc'^^^)) along the critical line 
P = Pi'^'-' to monitor the effect of quantum fluctuations 
enhanced by inhomogeneity. To this end, we perform a 
scaling study of the staggered magnetization on perco- 
lating clusters for fixed size Nc, with Nc up to 2048 and 
(3 = SNc- The obtained data are extrapolated to the 
thermodynamic limit by a three-parameter polynomial 
fit (mliNc)) = (mlioo)) + a N~^^^ + b N^. 

At the homogeneous bond percolation threshold Pd = 
Pl = 0.5 the percolating cluster is found to have anti- 
ferromagnetic long-range order, in agreement with Ref. 
0. However, the finite value of the order parameter is 
strongly reduced as the inhomogeneity is turned on, both 
in the "dimer" {0 ■k/2) and the "ladder" [9 ^ 0) di- 
rections, until it vanishes at non-trivial 6 values, 6l ~ 
0.34 in the ladder limit and 6'd ~ 1.23 in the dimer 
limit. This means that inhomogeneous percolation ex- 
hibits non-linear quantum fluctuations that are able to 
destroy the antiferromagnetic long-range order, in con- 
trast to the homogeneous percolation case-- with a renor- 
malized classical ground state. 

According to Eq. 10 this must be reflected in a pro- 
found change of the critical properties, since the van- 
ishing of the staggered magnetization on the percolating 
cluster decouples the magnetic transition from the perco- 
lation threshold. Beyond the two critical values 9l and 
6o, a higher concentration of bonds will be required for 
the system to magnetically order than for the lattice to 
percolate, thus opening up an intermediate phase with a 
novel quantum-disordered ground state. 
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FIG. 5: (color online) Upper panel: scaling of the disorder- 
averaged squared staggered magnetization on the percolating 
cluster with fixed size Nc- The different curves correspond 
to various points along the classical percolation transition. 
The continuous lines represent quadratic fits. Lower panel: 
extrapolated thermodynamic values of the staggered magne- 
tization from the upper panel as a function of the inhomo- 
geneity parameter 6. The dashed line is a guide to the eye. 



V. PHASE DIAGRAM 

In this section, QMC results are discussed for the mag- 
netic phase diagram in the {Pd,Pl) plane. The magnetic 
transitions are located via scaling of the spin-spin correla- 
tion length ^, extracted from the structure factor through 
the second moment estimator—. The inhomogeneity of 
bond dilution breaks the discrete rotation symmetry of 
the lattice, such that two distinct correlation lengths 
and £,y need to be considered along the x and y lattice 
directions. At critical points separating 2D long-range 
order from disordered phases, both correlation lengths 
must scale linearly with the lattice size, ^x,£,y ^ L, re- 
gardless of the universality class of the transition. In 
locating the magnetic transitions we have verified that 
this condition is satisfied. This leads to the conclusion 
that, despite its inhomogeneity, the magnetic system has 
unique well-defined phase transitions. 

The resulting phase diagram is shown in Fig. 0] The 
magnetic transition line coincides with the percolation 
line for moderate inhomogeneity. As discussed in the 
previous section, in this region of the phase diagram the 
percolating cluster at the percolation threshold is antifer- 
romagnetically ordered. Therefore not only the magnetic 
transition and the geometric transition coincide, but also 
the critical exponents of the magnetic transition are those 
of 2D percolation^^, according to Eq. ©. 

For strong inhomogeneity. the scenario of a classi- 
cal percolation transition in the magnetic Hamiltonian 



is precluded by non-linear quantum fluctuations. In fact, 
two multicritical points occur, in the ladder and in the 
dimer direction, beyond which the magnetic transition 
decouples from the percolation threshold. The location 
of these multicritical points is quantitatively consistent 
with the vanishing of antiferromagnetic order on the per- 
colating cluster at two critical values of the inhomogene- 
ity parameter 9 (Fig. [SJ. 




(a) (b) 




(c) id) 

FIG. 6; (color online) Percolating clusters for the bond- 
diluted square lattice with size 24 x 24 and periodic boundary 
conditions: (a) Pd = 0.15, Pl = 0.8; (b) Pd = 0.5, Pl = 0.5 
(homogeneous case); (c) Pd ~ 0.9, Pl = 0.45; (d) Pd ~ 1.0, 
Pl = 0.44 (randomly coupled dimers). 

Beyond the two multicritical points, intermediate 
quantum phases appear, in which antiferromagnetic or- 
der is absent on the percolating cluster both at and away 
from the percolation threshold. These regions represent 
novel quantum-disordered phases on an infinite family 
of 2D percolated random lattices. Representative struc- 
tures of this family are shown in Fig. |S1 together with 
the homogeneous percolating cluster. It is remarkable 
that the homogeneous percolating cluster (Fig. Elb)) 
has long-range antiferromagnetic order in the thermo- 
dynamic limit, whereas the other percolating clusters in 
Fig. EIa),(c),(d) do not. In fact, the latter structures 
have a significantly higher average number of bonds per 
site than the homogeneous percolating cluster, which has 
a fractal dimension^^i lower than 2 {Df — 91/48). A de- 
tailed discussion of the nature of the quantum-disordered 
phase will be given in Sec. IVIII 

The evolution of the critical exponents upon tuning 
inhomogeneity has been investigated in Ref. l23l Surpris- 
ingly, the two multicritical points do not mark an evi- 
dent discontinuity in the universality class of the model, 
and the classical percolation exponents appear to persist 
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also for intermediate inhomogeneity beyond the multi- 
critical points, changing to radically different values only 
in the extreme inhomogeneous limits. Therefore the mag- 
netic quantum phase transitions close to the multicritical 
points appear to retain a percolative character. This re- 
sult can be understood within the following argument. 
To destroy the long-range nature of a random network 
it is sufhcient to cut a few links on its backbone. This 
implies that the non-linear quantum fluctuations lead- 
ing to quantum-disorder on the percolating cluster need 
only have short-wavelength components, in contrast to 
conventional quantum phase transitions in translation- 
ally invariant lattices. Indeed, strongly fluctuating states 
appearing on links made of segmented ladders or deco- 
rated chains (see Sec. have a markedly local effect 
on magnetic correlations, i.e. they magnetically decouple 
the portions of the percolating cluster connected through 
that link. This effect is equivalent to geometrically re- 
moving that link, thus leaving two portions of the cluster 
disconnected^'^. If weak links occur on the backbone of 
the percolating cluster, this mechanism is sufficient to 
lead to quantum disorder. Only quantum fluctuations 
with wavelengths comparable to the length of the weak 
link are required. To reestablish magnetic correlations 
between the two portions of the percolating cluster dis- 
connected by quantum fluctuations, it is then necessary 
to add more bonds to the cluster in order to find an alter- 
native path for magnetic correlations to spread over the 
cluster. Since this process is of geometrical character it 
endows the quantum phase transition with a percolative 
nature. 

The above picture of a percolative quantum phase 
transition breaks down for strong inhomogeneity. In this 
limit, the building blocks of the percolating cluster are 
the strongly quantum fluctuating substructures depicted 
in Fig. 12 and hence the quantum phase transition is 
driven by the competition between the energy scale of 
the gap above the ground state of such substructures 
and the energy scale of the spatially random couplings 
between the substructures. In a real-space picture, the 
competition occurs between the finite correlation length 
of the substructures and the characteristic length set by 
the spacing between the spatially random couplings. This 
competition will be discussed in detail in the next section. 
The quantum phase transition is completely disconnected 
from classical percolation of the lattice, and it is reason- 
able to expect for it to be in a different universality class, 
as concluded in Ref. [2^ 



VI. CRITICALITY NEAR THE LADDER LIMIT 

So far, the onset of the quantum disordered phases 
was discussed, approaching them from the homogeneous 
limit of Eq. . A complementary understanding of the 
novel quantum phase and the associated quantum phase 
transitions can be gained by starting from the extreme in- 
homogeneous limit of randomly coupled ladders (9 = 0). 



This limit lends itself to a very simple analysis in terms 
of the properties of the single spin-1/2 antiferromagnetic 
Heisenberg ladder. 

In the limit of Pj, = 1, the system is composed of an 
array of two-leg ladders randomly coupled by D bonds. 
The probability of activation of a dimer bond is reflected 
in the average distance (r) of two D bonds along a column 
separating two neighboring ladders, namely Pd = ^1 {t)- 
A single D bond is sufficient to connect two ladders ge- 
ometrically. In other words, (r) = L is sufficient to con- 
nect all ladder subsystems on a L x L lattice. Taking the 
limit L — > oo, one obtains immediately Pd — *■ 0. Thus 
geometrically the system percolates at any infinitesimal 
concentration Pd > and the percolating cluster is the 
entire lattice itself. 




'ladder 



FIG. 7: (color online) Cartoon of the ordering mechanism 
in randomly coupled ladders. When Pd = « l/^iaddor 

(see text) percolating strings of block spins (rectangles) ap- 
pear in the system. The divergence of the correlation length 
along such strings drives the onset of long-range correlations 
between the strings and ultimately of 2D long-range order. 



However the system does not develop antiferromag- 
netic order until a finite value Pd ~ 0.32 is reached. 
This is not surprising since each two-leg ladder has a 
quantum-disordered ground state with a finite excitation 
gap. This gap has to be overcome before the Neel or- 
dered state can become the ground state. A more quan- 
titative argument can be formulated based on the fact 
that each isolated ladder has a finite correlation length^i 
6addci ~ 3.19. One can imagine a renormalization group 
transformation that creates effective block spins from all 
the spins within a correlation volume in such a way that 
all block spins are perfectly uncorrelated with each other 
(see Fig. [7||. When turning on the D bonds, it becomes 
evident that, to establish long-range correlations in the 
system, it is necessary to have a string of block spins per- 
colating from one side to the other of the system. Given 
that all block spins are independent within the same lad- 
der, this is only possible if, on average, each block spin is 
connected by a 13 bond to both its right and left neigh- 
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boring block spin. This requires an average spacing of 
the D bonds of (r) « bladder on each column, which in 
turn leads to the following estimate of the critical con- 
centration of dimer bonds 



Pr 



6 



ladd cr 



(7) 



Despite the simplicity of the argument, this estimate is 
surprisingly good: Pd,c ~ 0.313 to be compared with the 
QMC result Pd,c = 0.32(1). 




FIG. 8: Correlation length of the bond-diluted 128 x 2 ladder 
compared with the average ladder length (/). The data for 
gladder are taken at /3 = 2048; comparison with data at /? = 
4096 (not reported) shows no significant deviation. 



Interestingly, this argument can be extended to the 
case of randomly coupled weakly diluted ladders, namely 
for Pf, < 1. The correlation length of the ladder remains 
the dominant length scale of the problem, as long as it 
is much smaller than the new length scale introduced by 
the dilution 



(0 



1 



(8) 



corresponding to the average length of the ladder seg- 
ments after dilution. 

To test to which extent the above argument is appli- 
cable we performed SSE-QMC simulations on a single 
bond-diluted two-leg ladder to determine the evolution 
of the correlation length upon dilution. The results are 
shown in Fig. |S1 For low doping concentrations, we ob- 
serve that the average correlation length of the ladder in- 
creases moderately. For higher concentrations of missing 
bonds, the average ladder length {I) becomes compara- 
ble to the correlation length, which then crosses over to 
a decreasing behavior simply reflecting that of {I). 

We now focus on the increasing tendency for small di- 
lution. This behavior reflects various competing mecha- 
nisms. The dominant effects having the highest proba- 
bilities {0{Pl) and 0{Pl)) are: 



eff 



(a) 



03) 



(c) 



FIG. 9: (color online) Dominant effects of bond dilution on 
a two-leg ladder: (a) local antiferromagnetic modulation due 
to a missing rung bond; (b) effective couplings between spins 
missing their rung bond; (c) enhancement of the singlet com- 
ponent (indicated by an ellipse) of the two rung dimers adja- 
cent to a missing leg bond. 



(a) The removal of a rung bond results locally in two un- 
coupled S — 1/2 free moments which introduce a local 
staggered modulation of the adjacent spins'^*'^^ within a 
correlation volume ^ladder (Fig- 13a)). This leads to a 
local increase in the antiferromagnetic correlations along 
the legs; 

(b) If two rung bonds are removed sufficiently close to 
each other, the induced local antiferromagnetic modula- 
tions of the ladder ground state "lock in phase" , leading 
to an effective coupling between the free moments, expo- 
nentially decaying with the distance in the limit of low 
dopingSSiiSi, 



Jcff-^ (-ly-^ exp [- 



ladder I 



(9) 



These long-range interactions also lead to an increase 
of the correlation length in the system through an 
order-by- disorder mechanism (Fig. Etb)); 
(c) The removal of one leg bond, on the contrary, leads 
to a local enhancement of the rung-singlet component 
of the state of the two adjacent rungs (Fig. Efc)). This 
mechanism weakens the correlations along the two legs. 

Effect (c) has twice the probability compared to the 
effects (a) and (b) because there are twice as many leg 
bonds as rung bonds. The non-trivial competition be- 
tween correlation enhancement (a)-(b) and correlation 
suppression (c) has the combined effect that in bond- 
diluted ladders, the correlation length always remains fi- 
nite down to zero temperature. This is to be contrasted 
with the case of site-diluted ladders, in which the corre- 
lation length diverges algebraically with decreasing tem- 
perature, ^ ^ ^ reported in Ref. 0. (Obvi- 
ously, even in the site-diluted case the average correla- 
tion length is eventually upper-bounded by the average 
length of the ladder segments). It is interesting to note 
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that for the case of site dilution, correlation-suppressing 
effects are not present beside the simple fragmentation 
of the ladder. 
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FIG. 10: (color online) Phase diagram of the QHAF on the 
inhomogeneously bond-diluted square lattice (see also Fig. 0J 
close to the ladder limit and compared with the single-ladder 
prediction Pd,c = l/6addor(-PL) (see text). 



The absence of diverging correlations in bond-diluted 
ladders is crucial for the occurrence of an extended 
quantum-disordered phase. The finite correlation length 
of a single ladder for Pl ^ 1 allows us to repeat the RG 
argument sketched in Fig. dfor the case of bond-diluted 
ladders. One needs a finite and non- classical concentra- 
tion of dimer bonds Pd.c to drive the system into a Neel 
ordered phase, directly related to the quantum correla- 
tion length of the system bladder, as long as it is signifi- 
cantly lower than the classical length {I): 



PdAPl) 



1 



6addcr(^'L) 



(10) 



The above estimate is a theoretical prediction for the 
quantum-critical curve Pc = {Pd,c, Pl.c), which can be 
directly compared with the QMC results, shown in Fig. 
nn The range of quantitative validity of Eq. EH is sur- 
prisingly large, extending down to Pl ~ 0.875, which is 
interestingly also the value where the quantum correla- 
tion length bladder (-Pl) starts to cross over towards the 
decreasing classical behavior (Fig. |HJ|. At this point the 
ground state of the ladder segments is strongly altered by 
the presence of the dimer bonds, and hence the argument 
leading to Eq. ((Tn|l breaks down. 



VII. THE QUANTUM-DISORDERED PHASE: 
CORRELATIONS AND GRIFFITHS-MCCOY 
SINGULARITIES 

In this section, we discuss the nature of the quantum- 
disordered phase. Ground-state properties and thermo- 



dynamic observables are considered which reveal the na- 
ture of the low- lying excitation spectrum. A fundamental 
aspect of this phase, as it is typical for disordered sys- 
tems, are large variations in the local properties of the 
system, reflecting the local geometric structure defined 
by disorder. In the previous sections, we have discussed 
how the presence of segments of ladders or decorated 
chains leads to the local formation of RVB /dimer-singlet 
states. At the same time, bond dilution of ladders leads 
to the appearence of S* = 1/2 local degrees of freedom, 
effectively interacting across the strongly quantum fluc- 
tuating regions. The coexistence of such different subsets 
of spins within the same system raises the question of 
what the global nature of the ground state really is. 



A. Short-range correlations 
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FIG. 11: (color online) Temperature and finite-size scaling of 
the static structure factor S{Tr,n) for a representative point 
in the quantum-disordered regime {Pd = 0.15, Pl = 0.85). 
Dashed lines are a quartic polynomial fit in T = /3~^ [see 
Eq. lllH l to extrapolate to the (5 —* oo limit (when required). 
Inset: /3 — > oo extrapolated values of the static structure 
factor plotted vs. system size, clearly showing saturation in 
the thermodynamic limit. 



We first characterize the ground state in terms of global 
correlation properties. The absence of magnetic order in 
this phase has already been discussed in Sec. IIVI Fig. ITTl 
shows the staggered structure factor S'(7r, tt) for a repre- 
sentative point {Pd = 0.15, Pd ~ 0.85) in the quantum 
disordered phase on the "ladder" side of the phase dia- 
gram in Fig. 21 The finite-temperature data for S'(7r,7r), 
obtained using /3-doubling, are extrapolated to the limit 
T — > by accounting for power-law temperature correc- 
tions up to fourth order: 



4 



S'(7r, tt; T) = ^(Tr, 7r)^^„ + > -5- 



(11) 
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FIG. 12: (color online) Left panel: temperature and finite- 
size scaling of the correlation length along the two lattice di- 
rections {^x and for a representative point in the quantum- 
disordered regime {Pd = 0.15, Pl = 0.85). Dashed lines are 
a quartic polynomial fit in T = (3^^ to extrapolate to /9 — > cxd 
(when required). Right panel: /3 — > oo extrapolated values of 
the correlation length(s) vs. system size; the dashed lines are 
polynomial fits up to second order in L^^ on the points for 
L > 48 to compensate for the finite-size effects. 
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FIG. 13: (color online) Finite-size scaling of the T = static 
structure factor and correlation length(s) for a point in the 
quantum-disordered regime {Pd ~ 0.15, Pl = 0.9) closer 
to the ladder limit. The dashed lines associated with the 
correlation length data are quadratic fits S,{L) — £,{oo)+a/ L-\- 



where the a^'s are fitting coefficients. The extrapolated 
T — Q values of ^(Tr, tt) are shown in the inset of Fig. 
as a function of system size. They clearly display satu- 
ration for L — > oo, proving not only that the system is 
disordered, but also that it has finite-range correlations. 
This fact is further reflected in the correlation length 



£,x{y) along the two lattice dimensions, shown in Fig. 
A fitting procedure similar to Eq. Hll|l is used to elimi- 
nate polynomial finite-temperature correction. The £,t=o 
values so obtained (right panel of Fig. 1121) suggest con- 
vergence towards a finite value for L oo, although 
the correlation length is large for the particular point 
(Pd = 0.15, Pd = 0.85), and, even for the largest con- 
sidered size {L — 72), ^ does not display full saturation. 
A polynomial fit £,{L) = ^oc+a/ L+b/ L'^ for L > 48 yields 
saturation values £,x.oo = 9(1) and ^j,,oo = 12.5(5), result- 
ing in a considerable correlation volume of ^ 100 sites. 
It is important to keep in mind that the finite correlation 
length in the quantum disordered system varies strongly 
with the degree of inhomogeneity, and it can become as 
small as ^ ^ 3 in the limit of randomly coupled ladders. 
For instance. Fig. ^| shows the finite-size scaling of the 
static structure factor and correlation lengths for a point 
closer to the ladder limit, Pd = 0.15 and Pl = 0.9. For 
this point one observes excellent temperature saturation 
of the data up to L = 64 for (3 < 16384 without the need 
of any fitting procedure. Moreover, the zero-temperature 
data reported in Fig. 1131 show a much better saturating 
behavior within the system sizes considered, leading to 
the asymptotic values £,x,oo = 2.5(1) and ^y^oo = 5.4(2). 

A clear picture of the short-range nature of correlations 
also emerges from real-space images of 



the local staggered structure factor. 



S'(7r,7r;i) 



J '0' 



• the correlation function C(0,i) = {SQSf 
a reference site and all the other spins; 



(12) 



between 



• the energy of each bond -E(ij) 

all shown in Fig. 1 141 for a given sample corresponding to 
a representative point of the quantum-disordered phase 
on the "ladder side". Here the symbol (...)(, indicates the 
expectation value associated with the ground state of the 
specific sample considered, not to be confused with the 
disorder average (...). 

Starting with the bond energy, we observe that quan- 
tum fluctuations introduce a wide range of variability 
in its values, energetically promoting weakly correlated 
dimers or quadrumers and rung bonds on ladder seg- 
ments, at the expenses of the energy of the adjacent 
bonds (due to the fact that a. S = 1/2 particle can form a 
singlet with only one other S — 1/2 particle). In particu- 
lar the rare D bonds have in general the weakest energies, 
due to the strong tendency of the energy to be minimized 
on the ladder segments. This leads to a clear quantum 
suppression of correlations in the direction transverse to 
the ladders. 

The local staggered structure factor of Eq. l(T^ mea- 
sures the "effective correlation volume" (integral of the 
antiferromagnetic correlation function) around each spin, 
thus identifying the spins that are most strongly corre- 
lated with the remainder of the system. In Fig. I14r a') 
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FIG. 14: (color online) Real-space images of the bond energy 
(both panels) , local static staggered structure factor S{i) (see 
text) (a) and of the static staggered correlations C(0,i) be- 
tween a reference point on a dangling bond and the remain- 
der of the system (b), for the percolating cluster on a 32 x 32 
sample with Pd = 0.15 and Pl = 0.85 at ^ = 8192. The 
thickness of each bond is proportional to its energy, whereas 
the radius of the dots in (a) and (b) is proportional to the lo- 
cal value of S{i) and C(0, i). In (a) red boundaries highlight 
the regions where local correlations exceed the average value, 
namely where S{i) > 5'(7r, n). In (b) the red boundaries mark 
the correlation volume beyond which C(0, i) < 8 x 10"'', cor- 
responding to a distance (in units of the correlation length) 
|n|/C~ln[C(0,0)/C(0,i)l f^3.4 . 



we observe that only a portion of the spins has a sizable 
local structure factor, when they belong to regions with a 
larger local coordination number. In Fig. I14r b) we have 



picked a site on a dangling bond attached to a segment of 
a ladder, and calculated its correlations to the rest of the 
system. Compared with Fig. ICT a). we observe that the 
correlation volume around the chosen spin is considerably 
smaller than the cluster of "correlated" spins identified 
by the local structure factor. This implies that spins that 
are most correlated with the remainder of the system are 
not all correlated with each other. The finite-range na- 
ture of correlations is also clearly shown in the picture, 
along with the non-monotonic decay of correlations with 
distance. The latter is due to the fact that a dangling 
spin is mostly correlated with spins not involved in a local 
singlet state, even over long distances, due to the eflFec- 
tive long-range couplings discussed in Sec. IVII Nonethe- 
less, even in the presence of such couplings, correlations 
decay significantly with distance. This picture is to be 
contrasted with the case in which long-range couplings 
give rise to a genuine order-by-disorder phenomenon, as 
for instance in the site-diluted dimerized systems treated 
in Ref. \^ There the real-space image of C(0, z) shows 
very little spatial decay of the correlations between the 
free moments of the system. 



B. Griffiths-McCoy singularities 

In this section, the nature of the low-lying excitation 
spectrum in the quantum-disordered phase is discussed. 
In a clean system, the presence of a finite real-space cor- 
relation length implies finite correlations also in imagi- 
nary time, i.e. the existence of a gap in the excitation 
spectrum. In disordered systems, however, this is not 
necessarily the case. The disorder considered here is com- 
pletely uncorrelated in real space but perfectly correlated 
in imaginary time, such that long-range correlations in 
the time dimension (and thus absence of a gap) may co- 
exist with short-range real-space correlations. 

The nature of the low-lying excitation spectrum is re- 
fiected in the low-T uniform susceptibility. The global 
uniform susceptibility of the system is typically domi- 
nated by the Curie contribution from odd-numbered fi- 
nite clusters, which masks the behavior of the quantum- 
disordered percolating cluster. Hence, in the following 
discussion the contributions from the finite clusters are 
neglected. Instead, we focus on the uniform susceptibil- 
ity of the percolating cluster with an even number of sites 

Xu,c^^^Y.^StSI). (13) 

The results for Xu,c as a function of temperature, both 
in the ladder regime and in the dimer regime, are shown 
in Fig. 1151 A gapped spectrum would imply an expo- 
nentially vanishing susceptibility as T ^ 0. However, 
we find a non- vanishing susceptibility in this limit, and 
observe that Xu,c instead displays a power-law divergence 



11 



a=0.206(9) 
a=0.202(6) 
10 

a^O. 175(5) 
10' 

a=0.173(5) 




a=0.090(8) 
a=0. 110(5) 




• 




--0.85 


■ 




--0.9 


♦ 


Pd= 


--0.95 



lo- 



10" 



FIG. 15: (color online) Low-temperature uniform suscep- 
tibility of the percolating cluster in the quantum-disordered 
regime: ladder regime (upper panel) and dimer regime (lower 
panel). The solid lines are power-law fits of the form Xu,c ~ 
r~^+", and the resulting fit coefficients a are indicated. In 
the upper panel, the percolating cluster is picked as the largest 
cluster in a 64 X 64 lattice, whereas in the lower panel it is 
grown freely from a seed site up to a fixed site of Nc — 1024 
(see Sec. ITTTll. 



Xu.c 



(14) 



with a ^ 0, such that this divergence is not Curie-like. 
In particular, a — a{P) varies continuously as one scans 
across the quantum-disordered region: a « 0.09 0.2. 
This clearly shows that the system is not only gapless, 
but that its gapless excitations lead to singularities in the 
response to an external field. 

Quantum-disordered phases with divergent response 
properties at T ^ are familiar in the literature of 



disordered systems. They are classified as quantum 
Griffiths phases, and have been theoretically shown to 
exist in bond-disordered dimerized chains^, in bond- 
disordered S = 1 (Haldane) chains'*'^-'*^, and in disordered 
transverse-field Ising modcls"'^"'^^'^^ (see Ref . .45; for a re- 
cent review). As in our model, their central feature is the 
coexistence of short-range correlations together with so- 
called Griffiths-McCoy singularities in the response func- 
tions. Such singularities are due to the anomalous role 
played by local low-energy excitations living on rare re- 
gions of the system. 
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FIG. 16: (color online) Sketch of a low-energy excitation in 
the bond-diluted ladder system, (a) Four spins with miss- 
ing rung bonds are coupled across a clean region of length 
I with weak effective couplings Jos ~ exp[— o"?]; (b) rotating 
the spins at one end of the clean region by the same angle, 
spin excitations occur on the weaker bonds only, leaving the 
stronger ones unchanged. 



What are the excitations leading to the singular re- 
sponse in our model? To answer this question, we fo- 
cus on the ladder limit, referring to the discussion of 
bond dilution in ladder systems given in the previous 
section. There, it was observed that bond dilution leaves 
free S = 1/2 moments interacting via effective long-range 
couplings, Jeff, exponentially decaying with the distance. 
The presence of this weak energy scale in the system 
immediately populates the low-lying spectrum down to 
the ground state. For any arbitrarily small energy, one 
can always find two S = 1/2 moments which are suffi- 
ciently far from each other to be coupled with this en- 
ergy. Nonetheless, to be able to excite the system at this 
energy scale only, one needs the two moments to be sep- 
arated by a large clean region, such that one can rotate 
all moments at one end of this region of a given angle, 
leaving the higher-energy couplings at shorter range un- 
changed (see Fig. I16|) . Therefore on the one hand there is 
an exponentially small energy of the localized excitation. 
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FIG. 17: (color online) Probability distribution of the lo- 
cal susceptibility for a representative point in the quantum- 
disordered phase {Pd = 0.15, Pd = 0.9), taken at /? = 128L. 
The slope of the linear tail of the distribution is proportional 
to the a exponent in the thermodynamic limit (see text). 



leading to an equally weak local gap A: 
A exp[— cr/]. 



(15) 



but, at the same time, the presence of a large clean region 
of length I between the free moments is a rare event with 
exponentially small probability 



wil) ^ (1-Pl)' -exphc 



(16) 



By substituting Eq. I|15|) into Eq. 1)16(1 with the appropri- 
ate metric term to preserve the normalization condition 
of the probability distribution, the two exponentials com- 
pensate each other, leaving a power-law distribution^^ of 
the local gap A: 



u;(A) - A"/'"- 



(17) 



A very similar argument can be applied to the dimer 
hmit of the model. Here, e.g., an exponentially rare long 
clean decorated chain with weakly interacting spins at its 
ends leads to a similar probability distribution of local 
gaps. 

The above result has immediate consequences for the 
thermodynamics of the system^^, giving rise to the ob- 
served power-law behavior of the uniform susceptibility 
j^a-i -^yjjgj-g Q, = c/(T. In particular the constant c is 
clearly disorder-dependent [c — \ ln(l — Pi)| viaEq. ifTC)) ]. 
such that a is expected to be non-universal, as it is ob- 
served in our data. 

In order to verify the microscopic interpretation of the 
QMC results for Xu,c in terms of a power-law distribu- 
tion of the local gaps, Eq. (|17|l . we study an observable 



that directly probes the local gap, i.e. the local suscepti- 
bility^, 



'0 



(18) 



Hence, if the z-th spin is involved in a local quantum 
state with a local gap A, one obtains in the large-r (low- 
T) limit: 



(5f(T)5f(0))^exp[-Ar], 
such that, upon imaginary-time integration 



Xi, 



1 

A' 



(19) 



(20) 



By a simple change of variables in Eq. (|17|l one then finds 
a probability distribution for the local susceptibility, 

^(x.oJ^X-'-", (21) 
which is more conveniently recast in the form 

lnw(lnXioJ alnXioo: (22) 

i.e. the function lnw(lnXj^^) has a tail linear in (InXi^^) 
and a slope that directly yields the exponent a. 

Fig. 1171 shows results for the logarithmic distribu- 
tion of the local susceptibility, Eq. (|22|l . on sites of the 
percolating cluster only, for a representative point in 
the quantum-disordered phase {Pl — 0.9, Pd = 0.1) 
and for various system sizes at very low temperatures 
(/3 — 128L). Consistent with the Griffiths-singularity 
scenario, one observes a linear tail in In w(ln Xioc)' reflect- 
ing the power-law tail ofw{Xi„^) considered system 
sizes. This tail shows a strong size dependence, with a 
decreasing slope for increasing size, and crossing over into 
a non-linear finite-size tail for large Xioc values. The tail 
of the distribution w(Xioc) result of rare events as- 
sociated with clean regions and exponentially small gaps, 
such that finite-size corrections significantly affect the tail 
statistics. Nonetheless, the extracted values of the tail 
slope a converge towards the range a ~ 0.2 0.3, fully 
consistent with the independent estimate obtained from 
the low-temperature behavior of the global susceptibility 
Xu,c in the ladder limit. This result completes the picture 
of a quantum Grifhths phase for the quantum-disordered 
region of the system. 



VIII. CONCLUSIONS 

In this work, we have investigated the highly non- 
trivial interplay between quantum fluctuations and ge- 
ometric disorder realized in the inhomogeneously bond- 
diluted quantum Heisenberg antiferromagnet on the 
square lattice. We have observed that the inhomogeneous 
nature of bond dilution allows for a continuous tuning of 
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quantum fluctuations from linear to non-linear. Inho- 
mogeneity enables us to tune the system from a renor- 
malized classical phase, in which the magnetic transition 
is purely driven by geometric percolation, to a genuine 
quantum regime in which the magnetic transition is de- 
coupled from the percolation transition. Hence, a novel 
quantum-disordered phase appears between the magnetic 
and the percolation transition. This phase has the na- 
ture of a quantum Griffiths phase, characterized by finite- 
range correlations coexisting with a divergent uniform 
susceptibility. 

This represents, to our knowledge, the first evi- 
dence of a quantum Griffiths phase in disordered 2D 
Heisenberg antiferromagnets^, the previous observa- 
tions being limited so far to one-dimensional Heisenberg 
antiferromagnetgii2ii4 and one- and two-dimensional 
quantum Ising model3i2iiiii2ii&. Given the relevance of 
the two-dimensional Heisenberg model to the physics of 
cuprate oxides, we believe that this finding represents an 
important step towards the possibility of the experimen- 
tal observation of a quantum Griffiths phase. 

Obviously, the most natural form of disorder in real 
Heisenberg antiferromagnets is site dilution, realized 
by doping of non-magnetic impurities in the magnetic 
lattice^. Bond dilution can, in principle, be realized by 
vacancies in the non-magnetic lattice, in particular by 
vacancies on the sites of the non-magnetic ions involved 
in the superexchange paths. The presence of such va- 
cancies naturally alters the local electronic structure of 
the system, such that its effect on the magnetic Hamilto- 
nian may be more elaborate than simple bond dilution, 
and it might even introduce excess charge carriers in the 



system. The doping of the non-magnetic ions involved in 
the superexchange paths with other non-magnetic species 
introduces instead bond disorder—, whose effect can in 
principle be arbitrarily close to that of bond dilution. 

Beside the technical challenges in realizing bond dilu- 
tion (and bond disorder in general) in real antiferromag- 
nets, a fundamental aspect of this work is the central 
role of inhomogeneity in the search for novel quantum 
phases induced by disorder in two-dimensional Heisen- 
berg antiferromagnets. Recent studied show that the 
Neel ordered state of the square-lattice Heisenberg an- 
tiferromagnet is extremely stable towards homogeneous 
bond disorder, and homogeneous bond dilution destroys 
Neel order only at the classical percolation threshold'*^ 
Hence, inhomogeneity is an essential ingredient to realize 
strong enhancement of quantum effects through disorder 
beyond the previously investigated scenarios. 
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e.g., M. Continentino, Quantum scaling in Many-Body Sys- 
tems, World Scientific, Singapore, 2001). Nonetheless it is 
interesting to notice that the experimentally accessible ex- 
ponents, through e.g. neutron scattering, are the unshifted 
ones, (3 and v, along with the critical exponent of the static 
structure factor 7^ — — vz which concides with the per- 
colation exponent 7ci. 
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